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Abstract 

A scenario which overcomes the well-known cosmological overshoot problem associated with 
stabilizing moduli with steep potentials in string theory is proposed. Our proposal relies on the 
fact that moduli potentials are very steep and that generically their kinetic energy quickly becomes 
dominant. However, moduli kinetic energy red-shifts faster than other sources when the universe 
expands. So, if any additional sources are present, even in very small amounts, they will inevitably 
become dominant. We show that in this case cosmic friction allows the dissipation of the large 
amount of moduli kinetic energy that is required for the field to be able to find an extremely shallow 
minimum. We present the idea using analytic methods and verify with some numerical examples. 
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I. INTRODUCTION 



Let us start with a brief review of the framework in which the problem of stabilizing 
string moduli in the perturbative outer region of moduli space, where the string coupling is 
weak and/or the volume of the compact dimensions is large, is posed. 

Realistic models, for example in flux compactifications of string theory, usually include 
an effective N = 1 supergravity (SUGRA) theory below the string scale M s = 10 -2 (our 
conventions are such that M p = = 1.2 x W 18 GeV = 1) and supersymmetry (SUSY) 

breaking in some hidden sector at an intermediate scale Mj = 10~ 7 . General arguments 
based on symmetries show that the moduli superpotential must be a sum of exponentials in 
the moduli and perhaps an additional constant. These exponentials could be generated by 
stringy or field theoretic non-perturbative effects. 

In this framework, a typical potential for moduli fields a, has the N = 1 SUGRA form: 



Where K l i is the inverse metric derived from the Kahler potential K and Fj = diW + diKW 
is the Kahler derivative. Here we will consider a situation where all but one of the moduli 
have been stabilized at the string scale so that we can focus on the dynamics of one light 
modulus a as in the recent work of Kachru et al (KKLT) [| . 

The Kahler potential of moduli in the perturbative region is typically logarithmic in the 
fields, so that in terms of a canonically normalized component field the potential involves 
exponentials of exponentials. An important feature of such potentials is that they vanish or 
approach a finite constant when <fi — > 00. This corresponds to decompactification if <fi is the 
volume modulus and zero coupling if <p is the effective dilaton. This is the well known Dine- 
Seiberg problem Q. Therefore, if the potential has a minimum, it needs to be separated 
from the asymptotic region by a potential barrier. 

Realistic models need some fine tuning of parameters. SUSY breaking in the observable 
sector at the TeV scale requires that at the minimum of the potential F min = Mf = O(10" 14 ). 
In addition, if one assumes that the recent cosmological observations indeed indicate that the 
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and a typical superpotential has the form: 



W(a) = Y,Ae~ 



(2) 
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Figure 1: A typical moduli potential. The region of the shallow minimum had to be magnified 
by 26 orders of magnitude so that it can be seen. The vertical axis is in units of Mf, and the 
horizontal axis is in units of M p . 

cosmological constant (CC) is non- vanishing then at the minimum V m i n = O(10~ 120 ) > o\v\. 
so |F m j n | = v / 3|W / mm| + 0(1CT 12 °). A stable minimum further requires tuning of at least 4 
parameters to prevent tachyonic directions. Until recently it was very hard to find a single 
working model because the framework and parameters were too constrained. Now, with the 
development of models based on flux compactifications and the understanding of their vast 
parameter space, the discretuum, it has become possible to find models with minima in the 
outer region of moduli space jlj]. 

We are interested in estimating the height of the barrier that separates the minimum 
of the potential from the asymptotic region where the potential vanishes or approaches a 
constant. Let us assume that at some value <\> m in the potential has a true minimum. Since 
each exponential term is smaller in absolute value for <p > <p m ini we can generically expect 
that for > (fimin, \F\ < \F min \, and \W\ < \W min so that the height of the separating 
barrier is at most limited by the intermediate scale V max ~ Mf = 10~ 28 . If no further 
tuning is performed the height of the barrier is much lower than this estimate. Typical 
moduli potentials are therefore steep and have a very shallow minimum, when a stable one 
exists. A typical steep moduli potential with a shallow minimum is shown in Fig. d 

If one considers time-dependent solutions one encounters a cosmological version of the 
Dine-Seiberg problem, first discussed by Brustein and Steinhardt jjj. If the moduli start 




3 



from a generic point on the potential they are expected to reach the outer region of moduli 
space by classically rolling towards the asymptotic region. If they start to the left of the 
minimum, they will roll over the shallow barrier, and if they start to the right of the barrier, 
they will never reach the minimum and roll to the asymptotic outer region. To avoid this 
without additional sources, such as radiation, the initial position of the field, 4>q, has to 
be such that the initial height of the field, V((j)o), is at most an order of magnitude larger 
than the height of the barrier. Thus getting a bound solution requires fine-tuning of initial 
conditions to a very high accuracy. In addition, the steep potential inhibits the possibility 
of inflation while the moduli are rolling, since moduli kinetic energy tends to dominate the 
energy budget of the universe. The consideration of time-dependent solutions thus leads 
to additional criteria of stability of moduli beyond the standard static stability criteria, 
and therefore leads to additional requirements and constraints beyond the standard criteria 
which guarantee the absence of tachyonic directions. Such arguments have led us to propose 
previously that the most likely place for moduli stabilization is the central region of moduli 
space 

Several previous attempts to resolve the cosmological stability of moduli in general, and 
in particular the cosmological overshoot problem, have been made. Banks et al. have 
emphasized the role of the non-zero modes [5] and have noticed that they red-shift slower 
than zero modes. Dine |6 noticed the possible role of such slower redshift in stabilizing 
moduli. Barriero et al. j| discovered that the presence of additional sources helps to 



relax the problem, however not at the level that is required. Furthermore Huey et al. 
determined that if the moduli receive some specific temperature corrections then the other 
sources (ie. radiation) could be much more effective. We will comment on this further below. 
Of course all these discussions were made in the context of then known string models which 
were much more constrained than the models that we have now. 

Our proposed resolution relies on the existence of other sources, for example, a gas of 
relativistic particles: radiation. The nature of the additional sources and their energy density 
is not particularly important to us. A key fact that is crucial to the scenario that we propose 
is that kinetic energy (KE) in an expanding universe red-shifts at the fastest rate of all known 
sources. As the field rolls on its steep potential its KE builds up and very quickly leads to 
KE dominance. Since KE redshifts faster than the other sources, the additional sources 
will eventually become the dominant energy component. While these sources dominate, 
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they create a large amount of cosmic friction which dissipates a large amount of energy, 
allows the field to gently land into the basin of attraction of the shallow minimum, and to 
eventually settle down at the minimum. While these additional sources dominate, the field 
moves only a finite amount [3|>[3|- 

We do not rely on additional temperature dependent coupling beyond cosmic friction. An 
example of such coupling is T 2 2 . In fact, in the particular example that we use to illustrate 
our idea we assume that they are absent. If such couplings exist (as assumed by Huey et al. 
^]) they could further help in relaxing the constraints on moduli evolution. However, it has 
been argued that hig h temperature effects do not modify the form of the potential of string 



theoretic moduli 



Since our main goal here is to explain our idea and show that it can be realized rather 
than to explore in a general and systematic way the various possibilities and caveats, we 
concentrate on some explicit examples that allow us to discuss the basic argument. We will 
present a systematic search in a forthcoming publication. In Section II we use approximate 
analytic solutions to study modular cosmology with sources, and in Section III we verify nu- 
merically that our approximations and quantitative estimates are valid. Section IV contains 
our conclusions and a brief description of possible extensions of our scenario. 



II. MODULAR COSMOLOGY IN THE PRESENCE OF SOURCES 

We will discuss a specific model to expose the idea on which our proposal is based and 
to explain its basic ingredients. We consider a cosmology with a single field that has a 
potential of the form shown in Fig. ^ The field is assumed to start in the steep region of 
the potential. Some amount of radiation in a thermal state (with constant entropy so that 
T ~ I /a) is also assumed to be present p ra d = CT 4 = ^ where C and c are constants, T is 
the temperature and a is the Friedman-Robertson- Walker (FRW) scale factor. We further 
assume for simplicity a spatially flat universe. As we have explained, our proposal is not 
particularly sensitive to the nature of the additional sources. In fact all that is required is 
some additional source that redshifts slower than kinetic energy (see below). Here we will 
just incorporate radiation - in a later paper we will consider more general sources in detail. 
The qualitative features that we wish to illustrate will however remain the same. 

We do not discuss here the evolution prior to the epoch in which we can treat the effective 
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dynamics as a single scalar field in a cosmological background, or whether the universe starts 
in a quantum region. In both of these cases the universe would arrive at a starting point that 
is similar to the one that we assume. Such initial conditions can be arrived at in different 
ways. For example, a short period of inflation, a phase of so-called pre-big-bang evolution, 
or nucleation from nothing. 

The equations of motion that we need to solve are therefore 



Here a dot represents a derivative with respect to time, and H is the Hubble parameter 
H — a/a. 

Our scenario includes four distinct epochs that are described below. Each epoch starts 
with some specific initial values for the Hubble parameter and the field and its derivatives. 
The end values of the previous epoch provide initial values for the subsequent epoch. For 
definiteness, at the end of stage i we denote the values of the relevant variables by the 
subscript % so for example t<i is the time at which epoch number 2 ends. 

• epoch 1: Potential domination 

The field starts with zero initial velocity and a large potential energy on the steep 
part of the potential. The universe expands at a fast rate. If a substantial amount of 
radiation (or other sources) is initially present, the radiation energy density quickly 
redshifts as a -4 , and the main source of energy becomes the potential energy of the 



H 2 




(3) 



field. 



The equations of motion can be approximated by 



+ 



dV 

90 








(4) 



The energy E — ^- + V((p) is conserved in this epoch so all the potential energy is 
converted into the field's kinetic energy (KE). 
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The solution of the approximate equations (j3J) with the initial conditions at t — 0, 
4> = 4>o, V(<po) = Vq, and <j) = is the following, 



= y/2 [V - V(<f>)] 

t = 



p [Vq ~ V {<!>)] 



1 f 4> dcpJv^j)) 
In a = ' 



VQUo yJV -V {(/>) 

The velocity of the field at the end of this epoch 

0i = ^/2 [Vo-Vfa)], 



(5) 



(6) 



can be quite large if the potential is steep V ^> V(<pi), and the KE of the field becomes 
the dominant energy component. 

epoch 2: Kinetic energy domination 

The equations of motion in this epoch can be approximated by 



+ 3H(f) = 
1 „ 
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(7) 



In epoch 2 the equation of motion for <j) can be implicitly solved <fi = (j)\a\/ a 3 , implying 
that KE = |0 2 redshifts quickly KE = \^\a\/a & . 

The solution of the approximate equations (J7J) with the initial conditions at t = t\, 



i, and <j) = <pi = J2 [Vq — V(<j>i)], given in eq.(jHJ) is the following, 



In 



y> x {t - to + 1 



a(t) 



-a?0i(t - h) + al 



1/3 



(8) 



We can relate the amount of KE dissipation to the displacement in the field during 
this epoch 

(9) 
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which means that to dissipate a few tens of orders of magnitude in KE, as required, the 
field needs to move quite a bit in (reduced) Planck units to the extreme outer region 
of moduli space. If the field moves a several Planck lengths, as is the case in models 
that we have considered, then only a few orders of magnitude of KE are dissipated, 
approximately one order of magnitude per planck length displacement of the field. 

As the universe expands, both the KE and the radiation energy density redshift. Since 
the radiation redshifts at a slower pace, it will eventually "catch up" with the KE no 
matter how small its initial value. A possibility is that the potential energy will become 
dominant before the radiation. Whether this happens or not depends on the details 
of the model. In the successful cases (in the case of steep potentials) the radiation 
becomes dominant first. 

epoch 3: Radiation domination 

Epoch 2 leads to a radiation dominated epoch when the radiation "catches up" with 
the KE, (p vad ) 2 > (KE) 2 , 

~2 ±2 „6 



- » ^ (10) 



In epoch 3 the KE continues to redshift as K = and the expansion of the universe 
is faster, hence in this epoch KE is dissipated in a more efficient way. 

The equations of motion in epoch 3 can be approximated by 

4> + 3H<p = 



lc 2 



H = — A - (11) 



3 a 



The solution of eqs.lfTTj) is given by 



a = v v!^ - ^ 

= 02+ vgM_vgV! i (12) 

The displacement of the field during this epoch can be expressed after some algebra 
as follows: 



(, - ^ < V6. (13) 



\ (Prad) 2 
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Therefore, even if this epoch continues indefinitely the field will only move a finite 
distance P, Q]. KE on the other hand, continues to red-shift as 1/a 6 , so contrary to 
epoch 2, the field can dissipate a lot of energy while staying almost constant. 

• epoch 4: Potential domination 

Since the radiation energy density redshifts quickly, eventually, the potential energy 
density no matter how small, will come to dominate when V((p) 3> ^. In this epoch 
the equations of motion can again be approximated by 




At this point the solution of this system depends on the value of 0, and whether it 
ended up in the basin of attraction of a minimum, beyond the last minimum, or still 
at a high point. If the transition from KE dominance to RD occurs at a point that is 
close enough to the minimum (and to the left of the maximum), then the field will be 
trapped even in a shallow minimum since it has lost a huge amount of KE. 

The conclusion of this analysis is that as the field acquires a large amount of KE the latter 
seeds the elements of its own destruction by the lurking radiation. In the KE domination 
epoch the field does not move much per order of magnitude of dissipated KE, and in the 
radiation dominated epoch its motion is clearly bounded. Thus if the potential is steep 
enough, so that enough kinetic energy is acquired in the initial stage, and friction becomes 
important early on, and the difference in elevation between the starting point and the 
minimum is large enough, then with generic initial conditions the field would be bound. 

Obviously the detailed realization of this scenario depends on the values of the parameters 
in the potential and the nature of the specific sources. 

III. NUMERICAL EXAMPLES 

To verify that our approximations are meaningful, and to check their range of validity, 
we have made a series of numerical investigations. We report here only about some of them 
to illustrate some features of the scenario. 
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For the numerical investigation we have used the following potential for the volume 
modulus used in 



(15) 



where a = ev 2 / 3< ^. This potential was derived using the superpotential (coming from flux 
contributions and non-perturbative effects) 



W(<r) = W + Ae- acT , 



(16) 
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Figure 2: A bound case. Shown in the top panel are fractional densities of potential energy (blue), 
kinetic energy (red) and radiation energy (green) as a function of time. KE becomes dominant and 
then the radiation. Shown in the bottom panel is the evolution of the field as a function of time 
ending in the shallow minimum of the potential. 
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and the classical tree level Kahler potential. The last term in eq. (fT5|) comes from the 
contribution of an anti-Dbrane (Dbar brane). We stress that we use this particular potential 
merely for the sake of illustration of our mechanism. Similar results follow for potentials 
which incorporate alternatives to the Dbar brane term. 

All our numerical investigations were done in the conventions that M 2 = p \ = 1. We 

P 57TCrJV 

present our results in Figures El through To construct the numerical examples in these 

Q, Log 10 (V) 




Figure 3: A bound case. Shown in the top panel are the various energy densities (color coded as 
in Fig. 12} as a function of the scalar field position. Shown in the bottom panel is the potential 
and the starting position of the field. Note that the scale is logarithmic and that the difference in 
potential energy between the starting point and the minimum is about 30 orders of magnitude! 

figures we used the following values for our parameters: a = 0.1, A = 1.0, d = 3 x 10 -26 , 
Wq = —2.96 x 10~ 13 . For these parameters the potential has a a true minimum at = 7.06. 
At the minimum the value of the potential (the CC) is 6.35 x 10~ 42 and the barrier separating 
the minimum from the asymptotic region is located at <fi = 7.18. The height of the barrier 
is 6.28 x 10~ 34 . 

To illustrate the effect of radiation on the evolution of a moduli field we present examples 
of a bound and an unbound case. These two examples differ only in initial conditions. To 
create these two examples we set the initial conditions as follows. For the bound case: Vq = 
{Prad)o — 5.20 x 10~ 5 M 4 . The corresponding initial value of the field is <fro = 2.99, and the 
velocity of the field vanishes initially. For the unbound case: Vq = (p ra d)o — 3.78 x 10~ 2 M 4 . 
The fractional energy densities created by our bound solution are shown in Figures El and 
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Figure 4: An unbound case: various energy densities (color coded as in Fig. EJ) as a function of 
time. The only difference from the bound case is where the field ends up at the end of the radiation 
dominated epoch. In this case it lands to the right of the shallow minimum and continues to run 
on the potential in a "tracking" solution. 

El The first depicts the value of the field and the fractional energy densities as a function 
of time, while the second shows the same fractional energy densities and the potential as a 
function of the field. The four epochs mentioned previously in this paper are clearly visible 
in both representations. The fractional energy densities created by our unbound solution 
are shown in FigHJ Note that a significant fourth epoch is not present in this example. The 
field lands to the right side of the barrier, where a scaling solution quickly takes over. 

The addition of radiation to our model created an extra window, (2.99,3.95), of initial 
conditions which lead to bound solutions. The other interval, (7.04, 7.18), remains relatively 
unchanged with the addition of radiation. The two intervals of bound solutions can be seen 
in Fig El 

IV. DISCUSSION AND OUTLOOK 

We have proposed a scenario for resolving the cosmological moduli stabilization problem. 
We believe that we have identified the basic ingredients of the required solution. The 
application of the idea may vary according to the detailed models of string compactifications. 

Obviously this preliminary investigation needs to be followed by detailed and quantitative 
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Figure 5: Window of allowed initial conditions for bound solutions. Shown in the top panel are 
the various energy densities as a function of scalar field position (color coded as in Fig. for two 
bound solutions with two different initial conditions, depicted by solid and dashed curves. In the 
lower panel we show the potential as a function of the scalar field. The green dots on the potential 
are the end points of the two regions of initial conditions that lead to bound solutions. One of the 
regions is around the minimum, and the other is way up on the potential. 

exploration, including various sources, a systematic study with a wider range of parameters 
and potentials, and a quantitative analysis of stable potentials. In particular, the dependence 
of the range of bound initial conditions on the parameters of the potential such as the width 
and height of the minimum and barrier. 

Having found cosmologically stable models, it is clear that string models become less 
constrained. This provides a new perspective on the approach to string model building. 
Previously, because the theory was highly constrained, the hope was that only a single 
model, or a single class of models will satisfy the necessary conditions. Now the theory needs 
to be constrained as much as possible by data and phenomenological bottom-up constraints. 

As a final comment we wish to point out that constructing models that lead to slow roll 
inflation in this framework is rather easy to achieve. In models with stabilized moduli, one 
needs to find some that have a flat enough barrier as first discussed in the models of natural 
inflation Q. Sncfl a .node! in which the inflaton is the a*,o„ abated wflh the vohnne 

n 

modulus (j) was recently discussed in |14|. More models of inflation from a flat maximum of 



the potential of stabilized moduli were introduced in Q, and Q- 

Additional regions in 
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the discretuum may allow the construction of other types of inflationary models. 
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